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Abstract 

Most research and experiments carried out on the movement of particle! in air have 
been concerned with the total mast transport or flux of material from a surface It has 
been established thai the total rate of transport can be attributed to three model of motion 
which are not mutually exclusive; saltation, surface creep and suspension. Here, aalta- 
tion is the most important since It Is a prerequisite for the other two. The empirically 
determined saltation rate* for wind generated panicle movement are. In many cases, not 
in accord. Investigators have either modified the constants In Bagnold's saltation for- 
mula or developed complicated empirical equations to determine the transport rales. 
These equations are not likely to give further insight into the physics involved in the 
transport mechanisms. An alternative procedure It to examine the various components 
of the transport rate and determine the parameters on which it depends A novel flow- 
visualiiaiion/digital image analysts technique hat been developed to determine indivi- 
dual particle trajectories, from which saltation rates can be computed. This technique 
involves recording Individual particle motions, which have been illuminated with a 
strobed light source, on 35 mm high-speed film, digitizing the film and using a high- 
speed computer to track the panicle trajectories. The transport rates are determined sub- 
sequently from the particle-trajectory data, 



ulrodaKtioa 

In recent years, wind generated particle transport or wind erosion hat been receiv- 
ing increasing attention (Barndoroff Nielsen. 1985) from scientists and engineers 
interested in the effects of toil erosion, snow drift, sand movement and other related 
natural phenomena on local environmental, agricultural and economic conditions. For 
example, fugitive emissions are generated from open-area sources such at coal or iron 
ore stockpiles by the wind. These emissions are either transported downwind, affecting 
the air quality of near-by air-sheds or fall out locally, acting at a hindrance to a clean 
work place. Wind induced erosion can alto seriously affect the agricultural yield of a 
region, as was the case this year in the prairie provinces. In both of the above situations, 
valuable resources can be lost if appropriate control strategies are not adopted. In order 
to develop such control strategies, an understanding of the physics associated with wind 
erosion it required. 

Previous studies (Bagnold, 1941, Chepil. 1945, Gillette. 1978) have identified three 
primary modes of particle transport; surface creep, suspension and saltation. Surface 
creep pertains to particles greater than 1000 urn in diameter, which are loo heavy to 
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Since new estimate* of asthma prevalence in children in Canada 
nngi from 10 to 20», and possibly higher. it is clear that a 
substantial part of the Canadian population may be already 
suffering the ill effects of low-level exposure to acid ws*s, 
particularly sOi. 

(Supported by Ont . Ministry of Envt. I Health a welfare Canada) 
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become lirbonve but which cin roll or slide along die aurface. Thij transport mechanism 
accounts for about 10% of the total mass movement. Suspension pertains to panicles 
that generally have diameters less than 100 pun. These panicles remain airborne for long 
periods of time and are carried along at roughly the fluid stream velocity. The dominant 
mass movement (about 50% to 75% of the total) occurs by saltation, which pertains to 
particles having diameters between 100 to 1000 |lm. In this mode of transport, particles 
never go into complete suspension; rather, they remain close to the surface (usually 
within about 5 cm) forming a relatively thin 'saltation layer". Within this layer the parti- 
cles travel in a series of jumps or hops of about 10 cm in length. A fundamental finding 
of the aforementioned studies is that suspension of particles cannot occur without salta- 
tion as a prerequisite. 

Investigators have basically used a mass-flux technique in which the mass uplifted 
from a surface is collected by a vertical array of traps and subsequently analyzed in order 
to develop empirical relationships for the vertical and horizontal fluxes as a function of 
the surface shear velocity (u.) and other parameters. The first empirical relationship for 
the saltation rate (Q) was proposed by BagnoM (194 1) as 

Q-ul U> 

This has been modified and extended {White, 1979; Gillette. 1978) to include various 
parameters and is the basis for most calculations to determining the erosion rate associ- 
ated with a surface. Agricultural engineers have developed an empirical soil-erosion 
model, which involves a set of parameters, to predict the potential annual soil loss, 

E = f(I',C,K',L'.V). (2) 

where 1' is the soil erodibilily, C is the local climatic factor, K' is the soil ridge rough- 
ness factor, L* is the field length and V is the vegetation factor (Cbepil and Woodruff, 
1963). The above equation serves to illustrates the empiricism involved in determining 
the erosion rale. 

The mass-flux approach generally yields good engineering results for the quantity 
of material transported by the wind but does not provide any understanding of the phy- 
sics of the wind erosion process. An approach which does so is to consider the saltation 
process in terms of the motion of individual saluting particles within a turbulent 
boundary-layer in a Lagrangian frame of reference. Such an approach has been adopted 
by the Turbulence Research Group (TRO) in the University of Toronto's Mechanical 
Engineering (UTME) Department This involves the use of a novel flow- 
visualization/digital-image-analysis technique for obtaining the characteristics of indivi- 
dual panicle trajectories. In the present work, the technique used to obtain the trajec- 
tories is described, the governing equations of motion for the particles are discussed, and 
results for the particle flux profiles within the saltation layer, the saltation rate and the 
vertical mass flux of material are presented. 



Experimental Details 

The acquisition of good quality images or pictures is of crucial importance in the 
implementation of the aforementioned technique. The elements required for this include 
a wind tunnel facility, a lighting system, photographic equipment and film, and a digital 
image data acquisition system. 
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Experiments were carried out in a blow-through open circuit wind tunnel located in 
UTME, This tunnel has a rectangular test-section (180 cm long, 31 cm wide and 61 cm 
high) with a plexiglass top for transmitting light and a plexiglass side for viewing. 
Mono-dispersed (195 pm) silica sand particle beds. 100 cm long, 10 cm wide and 0.5 cm 
thick, were placed along the test-section centcrlinc. The entire test section was shrouded 
in black, except for a 25 cm by 25 cm viewing port located 40 cm from the leading edge 
or the panicle bed and a 10 cm slit directly shove the viewing port located along the 
centerline of the test section. A 35 mw He-Ne laser was used as a continuous eollimated 
light source. The laser beam was aligned above and along the centerline of the test- 
section and the beam deflected 90* by a prism and passed through a cylindrical lens (4.4 
cm focal length) to generate 1 12 cm by 0.2 cm plane of very intense light. 

Still photography was chosen in preference to high-speed cinematography, since 
with the former it was possible to capture most of the trajectory of an individual panicle 
in a single frame. As a consequence, it was easier to observe the particle's path, and the 
amount of digitization required was reduced. In order to resolve rapid changes in parti- 
cle position and velocity, the light source was "strobed" at 300 Hz by means of a rotating 
disk with an equal number of open and closed sectors. 

A 35 mm Contax 139 SLR camera fitted with a 105 mm telephoto lens mi posi- 
tioned perpendicular to the Row, outside of the test-section on a tripod. The distance 
font the plane of the film to the plane of light was approximately 41 cm, giving a field of 
view of 7.5 cm by 5 cm and a depth of view of 1 mm at an aperture opening of f2.8. 
Kodak Recording 1000 ASA film was used as the recording medium. The camera 
shutter was held open for 0.5 s and. with the particles illuminated by the stroboscopk 
light, "dashed" particle trajectories were recorded. These enabled measurements of veto- 
dry and position to be made in a Lagrangian frame of reference. 

A pilot tube was situated above the field of view to monitor the wind speed (5.5 
m/s) during etch experiment. Once particle motion was initiated, a roll of 36 frames, 
each containing a 0.5 s exposure, was shot every 2 seconds with the aid of a remote 
shutter cord and winder. Instantaneous wind velocity dan were acquired without the 
particle beds present using constant temperature cross-wire hot-wire anemometers. 
These data were recorded and subsequently computer-analyzed to obtain the mean longi- 
tudinal and vertical velocity components, turbulent uttetisities, velocity spectra and Rey- 
nolds stresses. 

After all the film was developed, the negatives were digitized on the UTME Digital 
Image Analysis System. This system has a Cohen CCD scanner with a Panasonic 25 mm 
T.V. lens, attached to a Quantex video processor, which is in turn connected to a tnicro- 
Vax computer. The scanner converts the image into 512 analogue scan lines which are 
digitized at a rate of 5 12 points per line by the Quantex, and the digital data are retrieved 
by the computer. Thus, each image Is represented as a 512 by 5 12 matrix with each pic- 
ture element or "pixel" containing 8 bits of tight intensity values, giving grey-levels 
between and 255. The digital data are subsequently stored on magnetic tape. 

Image Processing and Analysis 

A typical digitized image is presented in Figure I, where the dark dashed lines on 
the light background represent ibe panicle trajectories. Prior lo the extraction of the 
dashed- line data, the images had to be enhanced lo remove the noise (i.e., the uneven 
shades of background light and thin venical lines due to the plexiglass) while still 




151 



L 



1 



152 



maintaining die integrity of the image Due to the luge number of pinkies present and 
their various orientations, an * artificial intelligence/computer vision" approach wis 
needed 10 objectively analyse the pictorial data. Conventional approaches, such as thres- 
holding the grey levels or their gradients (Rosenfeld and Kak, 1982). did not provide 
adequate image enhancement, as is evident from Figure 2. Here, a modified gradient 
operator was used to enhance and delect the particle lines (Ctccone et al., 1987). 

A sophisticated scheme to further enhance the data was developed by considering 
the image as a two-dimensional signal and applying the Fourier transform to this signal. 
The dashed lines were found to have a particular "frequency" range which was different 
from that associated with the uneven background grey-level and spurious noise. A 
description of this scheme follows. 

The two-dimensional Fourier transform of an image, a complex function with real 
and imaginary components, is given by: 



F(u.v) = J / fdi.y^-^^'T'tfady, 



(3) 



where F is the Fourier transform of the image f, and u and v are the image frequencies 
associated with the x and y directions. Since the kernel function can be separated Into 
two pans, (3) can be carried out via two one-dimensional Fast Fourier Transforms (FFT) 
by first transforming the rows and then the columns of the image matrix or vice versa. 
Once the image has been transferred into frequency space, it can be readily filtered by 
means of a variety of digital filters. The "digital" representation of a band-pass filler is. 



H(u.v) = Htiu, v)H,Cu,v). 



(4) 



where H h (u. v) and ll,(u, v) are the high and low pass components of the filler. The 
frequency-domain form of the filtered image is given by, 
O(u.v) = F(u.v) H(u.v). 



This is transformed back into the real domain via an inverse Fourier transform, i.e., 



g(x.y) » J J G(u.v)e-' 2,,(u "**l'>dudv. 



(5) 



(6) 



In selecting the appropriate high and low filtering frequencies, an understanding of 
what contributes to each domain is required. Edges and other sharp transitions in grey 
levels are associated with the high-frequency -components while the uneven background 
levels contribute to the low frequency components. Thus, an image is sharpened by 
attenuating the low frequencies (high-pass filtering) and smoothed by attenuating the 
high frequencies (low-pass filtering). In the present study, a symmetrical Buoerwonh 
band pass filler was applied to the Image using 25 and 250 as the appropriate high and 
low cut-off frequencies. The result of applying the "Fourier" enhancement scheme is 
illustrated in Figure 3. As can be seen the integrity of the image is maintained, but the 
spurious noise has been removed and an even background of grey provided. It should be 
noted that this procedure is a robust and efficient means of enhancing an image. 

In order to locale the particle trajectories, which are thick lines, one could use a 
pattern-recognition technique to compare lines of various lengths and ingles with the 
enhanced image. This approach, however, would be extremely cumbersome and time 
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consuming. A more efficacious approach has been developed based on the fact that a 
line can be considered as a series of overlapping points, vis-a-vis the present situation . 
Thus, correlating the template representation of a solid circle with the image and subse- 
quently normalizing the resulting cross-correlation function produces a series of adjacent 
points which represent the centers of the dashed lines. The normalized cross-correlation 
function (Gonzalez and Wintz, 1977) Is given by: 



C,U,y)« 



fftto.P) g(x4<s.y+fl)dadB 
[/Jt I (a > B)dadB| , *[|fg 2 <x+«,y+f»dadp| v 



(7) 



where t represents the template and g represents an image. The template used ii shown 
in Table la. Since this template also responds to vertical lines, a second template of a 
verticil line (Table lb) was cross-correlated with the image and high responses to it 
were ignored. Figure 4 shows the local maxima of the above cross-correlation function 
and, as can be seen, the dashed lines have been extracted. The computer now knows the 
location of each pixel but does not know the plxel-to-pixel relationship. These pixels 
form a set of candidate points to be arranged into subsets which make up lines of finite 
length. Similarly, these subsets form the particle trajectories which must be linked 
together. This line tracking analysis can be considered a type of artificial 
intclligence/computer-vision analysis. 

Initially, the candidate pixels ate extracted from the digitized image column by 
column, i.e.. all candidates in column 1 followed by all ca n didates in column 2, etc. The 
following three tracking criteria are used: (i) the nearest neighbour (i.e.. Euclidean) dis- 
tance S V2; (ii) a minimum number of pixels is required to form a line; (iii) the pixels 
form a line which has no curvature. Criterion (i) is used to either initiate a line or to 
extend a line which is being tracked; (ii) It used to rid the image of residual noise; (iii) is 
used to break up lines which either cross or make sudden deviations, as is the case when 
a particle rises from or impacts upon a surface. This process is a combined parallel and 
serial approach where more than one line is followed at any one lime. Note that when 
lines cross, they arc flagged and any pixel which is accepted for one line is cross-checked 
with other line to see if it may belong to both lines. The tine subsets are then used to link 
the line segments together and reconstruct the trajectories. Figure S shows a typical 
result of applying the procedure. 

Since the light source was strobed at a 300 Hz, each panicle trajectory was divided 
into small straight line segments, with adjacent segments having similar characteristics, 
i.e., direction of flight and length. Thus, the difference in characteristics between dashed 
lines which belong to the same trajectory should be small. The line segments are linked 
together by first extending a given line by its own length in its direction of flight from 
left to right, i.e., the "end" coordinates of the extended line arc given by: 



x^-xj + Ax 



and 



(8.a) 
(8.b) 



y««i-yi+Ay. 

where x 2 and y 2 are the end coordinates of the original line of length AMAx'+Ay*) "* . 
Lines are linked together using the following three criteria: (i) the product of the slopes 
of any two lines must be positive or zero; (it) the difference between the start of a candi- 
date line and the end of the extended line must be less than a specified threshold; (iii) the 
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size of the lines must be similar. The lirst criienon is used lo ensure thai Hie lines are 
collinear and the second and third are used to ensure thai ihey belong lo the same panicle 
path. The result of using these criteria is presented in Figure 6. Note thai the panicle 
velocities are obtained by multiplying die length of the line segments by the strobing fre- 
quency. As illustrated, the panicle trajectories are incomplete. In order lo complete the 
trajectories, ■ Lagrangian random-flight model, based on the panicles' equations of 
motion is used. 

Lagrangian Random-Flight Model 

In the present case, the panicles can be considered spherical, and the requisite equa- 
tion of motion can be derived from a force balance on an Individual particle as follows, 
Newton's second law yields: 



dV. 



■EF.. 



m 



where m is the particle mass, V ii [he particle velocity and IF, is ihe sum of the 
relevant forces comprising drag, lilt, pressure, gravity and virtual mass as given below. 

mL (9.., 



F* M = -Qi7D 3 p r - 



FHa--QjD J p f (V»xa) 
k > <>V, 



■ -G>p-Pi)i 



F»i«Ml m*l * C„, — IT fJf 



dVl 
dl 



In equations (9.a) lo (9 e), C,, C, and C M are the coefficients of drag, lift and virtual 
mass, respectively, g it gravity, D is the panicle diameter, p and p ( are the panicle and 
fluid densities and V R is the panicle relative velocity given by: 



V.-fV,, 



+ V,, 1 )". 



Vu-V„ 



-u. 

-V, 



where U and V arc Ihe components of Ihe fluid velocity (V,) in the x and y directions, 
respectively. 
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After substitution and rearrangement, Ihe resulting equaiion is 

1 ^l 



Since p » p. ihe above equaiion can be reduced to the following, with only the 
drag, lift and gravity forces being retained: 



(12) 



The drag 



(13) 



dl 4 p, D p, 

This equaiioa provides a Lagrangian description of ihe panicle path, 
coefficient C f varies with the particle Reynolds number 

|V«|D 

and it calculated from empirical relationships (Morsi and Alexander, 1972). The lift 
force Is given in icnra, of the vector product of the relative panicle velocity and the fluid 
vorticily, 

(i»VxV,. (M) 

and was considered by Hunt et al. (1984). This force is similar to (he lift force which 
acts on an airfoil in an invfscid Bow with circulation. 

Given an initial panicle velocity and position and the wind field, (II) can be 
integrated numerically, both forward and backward la time, to determine the particle 
impact and ejection points and velocities. By randomly varying the wind field com- 
ponents and ensemble averaging the trajectories, a stochastic description of the 
panicle'! trajectory can be realised. Using the wind field measured by Ihe hot-wire 
anemometers, the trajectories are completed at shown in Figure 7. It should be noted, 
thai the fluid hat been assumed lo be unaffected by the presence* of the particles. 

Approximately 100 trajectories were recorded and analyzed to determine various 
physical characteristics such at the particle ejection velocity vector, maximum height 
and jump length. These trajectories arc consistent wilh those described by Owen (1964), 
While and Schuli (1977) and Hunt and Nalpanis (1985). A synopsis of the 193 p-m par- 
ticle retullt for a J J m/i wind are presented In Table 2. From the latter, It la apparent 
that the mean ejection velocities are greater than both u. and the fluid threshold velocity 
(u. ). This would indicate a mechanism other than mean shear stress acta so eject parti- 
cles in lo the boundary- layer. The most likely mechanism it impaction from other parti- 
cles rather than the lift force generated by the wind shear, rince thit force corresponds to 
less than 6% of gravity. 
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Estimation of the Mass Movement 

With the aid of the actual trajectories, estimates of the panicle fluxes, concentra- 
tions and saltation rates can be made, i.e. Lagrangian data can be used to obtain Eulerian 
point data. On the assumption that saltation is uniform and stationary, the mass flux uf 
particles across unit width of a plane perpendicular to the fluid stream is 

i 

G(D)=fmxV 10 dy. (15) 

where G is the saltation rate, x is the particle concentration, m is the mass of a panicle of 
diameter D, V is the particle horizontal velocity and y* is the height of the trajectory 
traced out by the particle (Owen, 1964). The integrand represents the horizontal flu* of 
material Now, (IS) can be rewritten as. 



0(D) = m«)(D)U(l;D)dt 



(16) 



where *(D) is the vertical flux of panicles, i(t;D) ii the particle horizontal velocity as a 
function of lime and I* is the time of travel. This results in: 

G(D) = mt(D)x"(D), (17) 

where x(D) is the particle Jump length (Jensen and Sorensen, 1982). 

Evaluation of (13) and (17) using the calculated panicle trajectories is carried out 
by recording the number of panicles caught In an ideal "numerical trap" and their hor- 
izontal velocities. This trip is unit width, Ay units high and is located at the maximum 
height of a particle trajectory. The horizontal flux of panicles is determined at each trap 
height, and from this, the saltation rale is calculated using (IS). Also, since the mean 
panicle jump length is known, the vertical flu* is calculated from (17). 

In order to obtain stable concentration and flux profiles, the ejection velocity com- 
ponents of each panicle are varied several limes, using a Gaussian random number gen 
crator. particle jump lengths are determined and panicle data are ensemble-averaged 
until convergence is attained. It was found that repeating the procedure (I.e., varying the 
velocity components) 500 times provided adequate convergence. Figure 8 presents the 
normalized horizontal fluxes as a function of height obtained for the 193 urn panicles. 
The saltation rate was calculated to be 0.124 g/cm/s using (15) and the vertical suspen- 
sion rate was calculated to be 0.017 g/cnAs using (17). It Is noted that the trend or the 
present results is similar to that of the data presented by Hunt and Nalpanis (1985). 

Conclusions 

Flow-Visualization of saluting silica-sand panicles in a wind-generated boundary- 
layer has been achieved. The resulting panicle-trajectory data have been objectively 
analyzed using an artificial intelligcnce/computer-vision approach A two-dimensional 
Lagrangian model has been developed and implemented to complete the panicle trajec- 
tories and from this, the mass flux of material has been estimated. 

The panicle statistics, indicate that impacting particles have a substantial effect on 
the ejection of other panicles since the particle ejection velocities are greater than the 
surface friction velocity. But, as suggested by Owen (1964), the introduction of panicles 




uiiu the fluid boundary-layer results in a slowing down of the fluid velocity, since the 
panicles act as momentum sinks; thus, panicles that are subsequently ejected do not 
gain as much momentum. It is suggested that the fluid boundary-layer acts as a feed 
back mechanism. 

A hypothetical scenario for the saltation process would be as follows. 
(i) Initial panicles begin to saltaie when the wind exceeds the particles fluid threshold 

velocity. 

(ii) Pinkies gain momentum from the boundary layer and impact onto the surface; the 
boundary-layer adjusts itself due to this exchange of momentum 

(iii) Impacting panicles transfer momentum to the surface, in addition to that transferred 
via the boundary-layer. 

(iv) The momentum transmitted to the surface is partitioned between the impacting par- 
ticles which continue to salute and other panicles causing them to eject. 

(v) Steps (ii) lo (iv) are repeated until a steady-state condition is reached. 
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Table 1 J Template representation of a solid circle used 
in template matching with image. 
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Table Lb Template representation of a vertical line used 
in template matching with image. 
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Table 2. Mean trajectroy characteristics for 195 pm panicles 
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Figure I Digitized negative of 195 um panicles moving left to right 
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Figure 2 Sobel gradient operator applied to image. 
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Figure 3 Fourier enhancement of raw image wing a Bwierwonh band pas, filler (23 lo 230). 
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Figure 4 Image after applying template matching to enhanced image. 
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l-'igure 5 Results after applying line tracking procedure. 
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Figure 6 Resulis afier applying line tracing rouline. 
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